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ABSTRACT 



The morphology of the infrared (IR) emission and the properties of discrete far- 
infrared (FIR) sources in M31 disk are studied using the high resolution (HiRes) maps 
from IRAS. Very thin and bright FIR arm segments are shown in these maps, which 
have similar structure as the HI gas but with much enhanced arm/inter-arm contrast, 
typically a factor of ~ 5 on the 60/zra image. We identify 39 unconfused sources 
(excluding the nucleus) and 14 confused sources (7 pairs) at 60 /xm by direct Gaussian 
fittings to the image. IRAS colors of 10 bright isolated sources are studied, which are 
presumably representative of the discrete sources in general. The colors follow the well 
known anticorrelation between f60/i/fioo/i and fi2/i/f25/i- All sources coincide with optical 
HII regions. A comparison with Ha observations shows that the total luminosity (in the 
wavelength range 8 — 1000/zra) associated with HII regions is 7.2 ±2.9 1O 8 -L , namely 
30 ± 14% of the total IR emission of M31. Half of the IR luminosity of this component 
emerges from the sources. This luminosity translates into a present-day star formation 
rate of 0.36±0.14 Mq /yr, about an order of magnitude lower than that of the Milky 
Way Galaxy. 

Subject headings: galaxies: individual - galaxies: interstellar matter - galaxies: 
photometry - interstellar: grains - interstellar: HII regions 



1. Introduction 

Two characteristics make M31, the Andromeda Galaxy, unique for a far-infrared 
(FIR) study. The first is its proximity: being the nearest spiral galaxy outside the Milky 
Way, it offers the best linear resolution for observations of spirals with a given angular 
resolution. This is a significant advantage for FIR observations, given the relatively poor 
angular resolution of the current generation of FIR instruments. Second, it is known to 
have a very low current star formation rate, perhaps an order of magnitude lower than 
that of the Milky Way (Walterbos 1988). Because of their very low FIR luminosities, 
not many such quiescent spirals can be studied in the FIR, and none except M31 with 
such high spatial resolution. 

The FIR emission of M31 has been studied by Habing et al. (1984), Soifer et al. 
(1986), and very extensively by Walterbos and Schwering (1987, hereafter WS87) using 
IRAS maps. The main results from these early studies can be summarized as follows: 

(1) M31 emits only a few percent of its total luminosity at wavelength longer than 
10/zra. The rest is mainly the optical radiation from stars (Habing et al. 1984; 
WS87). 

(2) In the central part of the IRAS maps, corresponding to the optical bulge, the 12 
and 25/zra radiation is due to circumstellar dust emission from late type stars, 
while the 60 and 100 \im emission is due to the interstellar dust heated by the 
interstellar radiation field (ISRF) powered by the old bulge stars (Soifer et al. 
1986). 

(3) The disk IR emission is strongly concentrated in a ring with a semimajor axis 
of ~ 45' and coincides in detail with the distribution of HII regions (Habing et 
al. 1984; see also Devereux et al. 1994). When it is decomposed into a warm 
(40K) component associated with star formation regions and a cool component 
associated with 'cirrus' (Cox et al. 1986; Helou 1986), less than 10 percent is due 
to the former (WS87). 



1 



(4) The strong mid-infrared (12 and 25 /zra) excess and the rather constant 60-to-100 
\im ratio in the disk show evidence of the existence of very small grains/large 
molecules (WS87; Helou 1986; Xu & Helou 1994). 

Recently, Devereux et al. (1994) made a comparison between the FIR emission and the 
Ha emission of M31, and argued that the close correspondence between the FIR and 
the Ha images suggests a common origin for the two emissions. According to these 
authors, up to 70% of the FIR emission of M31 could be associated with massive star 
formation regions, directly conflicting to the earlier result of WS87 (see above item 3). 

We are carrying out a new FIR study on M31, utilizing the high resolution maps 
(of ~ 1' resolution) produced with the HiRes software (Aumann, Fowler &: Melnyk 1990; 
Rice 1993). Our motivation is to study the properties of FIR sources and the diffuse 
emission in the M31 disk separately. This was not possible in the earlier studies because 
of the relatively coarse resolution of IRAS maps. In an earlier paper (Xu &: Helou 1994), 
we studied the IRAS color-color diagrams of a complete sample of small areas (2' X 2') 
in the M31 disk, and found evidence for a deficiency of Very Small Grains (but not 
the Polycyclic Aromatic Hydrocarbon molecules) in the areas dominated by the diffuse 
emission. 

In this and a companion paper, we present a study on the general properties of FIR 
sources and of the diffuse emission in M31. In this paper, Paper I, we study the overall 
morphology of the FIR emission and the properties of the discrete sources in the M31 
disk. In order to address the controversy as to the extent to which massive ionizing stars 
are responsible for the FIR emission of M31 (WS87; Devereux et al. 1994), we estimate 
the ratio between the FIR and Ha fluxes of HII regions using those 60/zra source which 
are fully covered in the Ha survey of Walterbos and Braun (1992, hereafter WB92). We 
present the new HiRes maps of M31 in Section 2. The procedure of source extraction is 
described in Section 3. The FIR colors of some 'isolated' sources are studied in Section 
4. In Section 5 we calculate the total dust emission of M31, from the sources and from 
the diffuse component respectively, and an estimate of the present-day star formation 
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rate of M31 is made. Section 6 contains the summary. Throughout this paper, we 
assume for M31 a distance of 690 kpc (1' = 200 pc along the major axis), an inclination 
angle of i = 77° [i = for seen face-on), and P. A. = 37°. 



2. HiRes IRAS maps 



The high resolution IRAS maps of M31 at 12/zra, 25/zra, 60/zra and 100/zra are 
new IRAS products from the Infrared Processing and Analysis Center (IPAC), using the 
IPAC High- Resolution (HiRes) processor (Aumann, Fowler & Melnyk 1990; Rice 1993). 
After the 20th iteration of the deconvolution, the resolution achieved is 0'.9 X 0'.5 (in 
cross-scan and in-scan directions respectively) for the 12 and 25/zra maps, l'.O X 0'.8 
for the 60/zra map, and 1'.5 X 1'.5 for the 100/zra map (see also Rice 1993). However, 
the resolution is not uniform across the maps (Fowler &: Aumann 1993), and depends 
in particular on the surface brightness of the background. In order to overcome this 
problem, and also to simplify the comparison between the four maps, we smooth all of 
them to a 1'.7 circular beam on a grid with 0.5' pixels. The details of the production 
of the maps, and a comparison with previous M31 IRAS maps, can be found in Xu &: 
Helou (1994; see also Xu & Helou 1995, hereafter ERRATUM). 

These smoothed maps are presented in Figures 1-4. Allowing for the difference 
in sensitivity, the four maps have similar morphology, characterized by a bright nucleus- 
plus-ring structure with the ring coinciding with that of the HII regions (Pellet et al. 
1978) and the nucleus coinciding with the optical bulge. This has been noticed in 
previous studies (Habing et al. 1984; WS87). The new HiRes maps however reveal 
features with a linear scale of a few hundred pc to ~ 1 kpc, namely the bright sources 
(giant star formation regions) and the very narrow widths of bright arms. Indeed it 
is rather surprising that the arms are so thin that they are hardly resolved by the 
HiRes beams in the cross-arm direction even on the original HiRes maps (resolution 
of ~ 1'). Compared to previous results (e.g. WS87), the contrast between the arm 
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segments (including the 'ring') and the inter-arm region is significantly enhanced in 
the new HiRes maps (e.g. the 60/zra map), with a typical value of about a factor of 
5 (avoiding the bright sources). Unfortunately, one cannot exclude the possibility that 
the contrast is artificially exaggerated by HiRes. These results will however be tested 
in the upcoming ISO mission. 

In the four panels of Fig. 5, we compare the surface brightness distribution of the 
60/zra map and that of a 1.'7 resolution HI map (Brinks 1984; Brinks &: Shane 1984) 
along two cuts, running through the center of M31 at position angles of 24° and of 50°, 
respectively. It appears that on the two maps there is good correspondence between the 
positions of features (dips and bumps), and the bumps representing the arm segments 
have about the same widths. On the other hand the HI emission decreases toward the 
center of the galaxy, while the 60/zra emission does not. In comparison with optical maps 
of similar angular resolution (Walterbos &: Kennicutt 1987), the FIR nucleus is smaller 
than the optical bulge, and the disk emission is much less smooth in the FIR than in 
the optical. An anti- correlation between the optical and the FIR surface brightness can 
be found in some places, e.g. the dust lanes which appear as bright thin arms on the 
FIR map, a natural consequence of the dust absorbing optical radiation and re-emitting 
in the FIR. Devereux et al. (1994) compared the HiRes FIR (40 - 120/zra) map with 
the Ha CCD image obtained with the Case Western Burrell Schmidt telescope at Kitt 
Peak, and found "striking correspondence" between the two. 

3. Source extraction 

There are many bright discrete sources on the HiRes maps. The brightest ones are 
concentrated on the famous 'ring' at about 9 kpc (45') galactocentric radius. Compared 
to the HII regions map (Pellet et al. 1978; WB92; Devereux et al. 1994), these sources, 
except the nucleus, coincide exclusively with bright HII region complexes. A comparison 
between the 60/zra sources and the optical HII regions (WB92) reveals that individual 
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FIR sources include typically 3 to 5 bright HII regions. The companion elliptical galaxy 
M32 does not show up in any of the four IRAS maps (Habing at al. 1984). Bright 
background radio sources (Berkhuijsen et al. 1983), which do not belong to the M31 
disk, are not detected, either. 

Because of the coarse resolution of IRAS, no analysis on discrete FIR sources in 
M31 could be carried out in previous studies. Even with the HiRes maps, many sources 
are still affected by confusion, in particular on the 100/zra map where they are more 
extended (Fig. 4) and the spatial resolution is the poorest. The 12/zra and 25/zra maps 
have low signal-to-noise levels. Hence an overall source extraction is done for the 60/zra 
map only. We at first run an automatic source subtraction program which fits individual 
single elliptical Gaussians over areas of 4'. 5 X 4'. 5. For a source extraction threshold at 
a surface brightness of 1 MJy/sr above the background set by the local diffuse emission, 
we found 70 sources (excluding the nucleus) on the 60/zra map (Xu and Helou 1993). 
The limit was assigned because the one-cr random fluctuation of the diffuse emission in 
the 60/zra disk of M31 is of the order of 0.5 MJy/sr. Therefore the sources are reliable 
at ^ 95 percent confidence level (> 2a level). Because of the source confusion problem, 
the parameters (e.g. the major axis and minor axis, the flux, and even the detection) for 
many sources vary significantly when the size of the fitting area is changed. In order to 
assess the uncertainty due to this factor, we repeat the source extraction program using 
two more sizes of fitting area, i.e. 3' X 3' and 6' X 6'. Out of the 70 sources detected in 
the first run (with the size of the fitting area of 4'. 5 X 4'. 5), 39 are detected in both of 
these two latter runs. They are represented by crosses in Fig. 3. 

For each of the 39 sources we have estimated, and reported in Table 1, the values 
and uncertainties of the following parameters: the coordinates R.A. (1950) and DEC 
(1950), the peak value of the 60/zra surface brightness S p , the major (maj) and minor 
(min) axes, and the position angle of the major axis P.A. (measured from the North 
towards the East). For any parameter x, the value is estimated from the weighted mean 
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of the logarithms of corresponding results in the 3 different fittings: 



E Ui ln(xi) 

i=l 



< ln(x) >= ^ (1) 

E Vi 

i=l 

where X{ represents the value of the parameter found in one of the fittings, and Ui the 
weight : 

(o-i/xi) 2 

where the a{ is the uncertainty on Xj. The variance of the logarithm is estimated with 
the same weight, and the systematic uncertainties due to the different Gaussian fittings 
are allowed for by introducing a scaling parameter s (Eadie et al. 1977): 
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and 

™ n 



E ^i(^i/xi) 2 E <"i [ln(xi) - < ln(x) >] 



- 2 (ln(x)) = s x ^ = -^-^ = , (4) 

E Vi E Vi 

i=l i=l 

with n = 3. The values and uncertainties listed in Table 1 are calculated from these 
means and variances: 

3 <ln(x)> 



x = e <ln W> (5) 



<r(x) = x^ 2 (ln(x)) . (6) 

The uncertainties of the peak surface brightness and of the total flux of each source 
include also the calibration uncertainty of the 60/zra map, which is at 4% level (Xu &: 
Helou 1994). 

The columns of Table 1 are arranged as follows: 
Col.(l): source ID in the list; 

Col. (2): Line 1: the 1950 equinox right ascension (R.A.) in hours, minutes and seconds; 

Line 2: uncertainty on R.A. in seconds of time; 
Col. (3): Line 1: the 1950 equinox declination (DEC) in degrees, arcminutes and arcseconds; 



Line 2: uncertainty on DEC in arcseconds; 
Col. (4): Line 1: peak 60/zra surface brightness S p in MJy/sr; 

Line 2: uncertainty on S p ; 
Col. (5): Line 1: major axis of the source in arcminutes; 

Line 2: uncertainty on the major axis; 
Col. (6): Line 1: minor axis of the source in arcminutes; 

Line 2: uncertainty on the minor axis; 
Col. (7): Line 1: position angle of the major axis in degree; 

Line 2: uncertainty on the position angle; 
Col. (8): Line 1: 60/zra flux in Jy; 

Line 2: uncertainty on the 60/zra flux; 
Col. (9): Line 1: sum of the Ha fluxes of the corresponding HII regions in 
10 _12 erg/cm 2 /s, taken from WB92; 

Line 2: uncertainty of the Ha flux; 
Col. (10): ID's of the corresponding HII regions in the catalog of WB92 for ionized nebulae 

in the NE half of M31. 



In addition, another 14 confused sources form 7 pairs (Table 2). For each pair two 
sources are detected individually in the test run with the fitting area size of 3' as well 
as in the original run with the fitting area size of 4'. 5, but in the run with the fitting 
size of 6' the two components are not distinguished, and therefore only one large source 
is detected. These confused sources are marked by open diamonds on the 60/zra map 
(Fig. 3) and listed in Table 2. The combined flux ieo^ and its uncertainty is calculated for 
each pair using the similar formulae as given in Eq's (l)-(6) with $60/i,i = $60/1,1 +^60/i,i 
and $60/1,2 = $60/1,2 + $60/i,2i where subscripts 1 and 2 denote the values found in the 
Gaussian fittings with areas of 3' X 3' and 4'. 5 X 4'. 5 respectively, and superscripts a 
and b the values for the component a and b. The flux $60/i,3 is taken from that of the 
single source (including both components) found in Gaussian fitting to the 6' X 6' area. 
The ieo/i of source pair 4a+4b has a large uncertainty (78%) due to the large dispersion 
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TABLE 2. 

List of Confused 60/zra Sources (Source Pairs) 



(1) 


(2) 


(3) 


(4) 


(5) 


ID 


RA (1950) 


DEC (1950) 




unc 




h m s 


o / // 


Jy 


Jy 


la 
lb 

la+lb 


00 38 13.0 
00 38 21.4 


40 20 15.2 
40 21 47.6 


4.43 


0.34 


2a 
2b 

2a+2b 


00 37 56.7 
00 38 09.4 


40 27 04.4 
40 28 49.0 


2.89 


1.32 


3a 
3b 

3a+3b 


00 38 40.8 
00 38 50.1 


40 33 42.8 
40 35 10.0 


5.28 


0.79 


4a 
4b 

4a+4b 


00 40 00.8 
00 39 54.8 


40 43 47.4 
40 45 50.9 


2.76 


2.15 


5a 
5b 

5a+5b 


00 39 22.0 
00 39 09.4 


41 07 09.9 
41 08 16.2 


2.51 


1.13 


6a 
6b 

6a+6b 


00 41 60.0 
00 42 12.4 


41 11 07.1 
41 14 42.3 


9.27 


1.92 


7a 
7b 

7a+7b 


00 40 16.9 
00 40 24.9 


41 20 53.6 
41 21 57.5 


9.65 


0.90 



among the 3 fittings (the factor s Eq(3) is large). Because of the confusion problem, 
for individual components of each pair we list only the right ascension and declination 
which are calculated from the weighted mean of the corresponding values found in the 
two Gaussian fittings using areas of 3' X 3' and 4'. 5 X 4'. 5 respectively. 
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4. IRAS colors of isolated sources 



From Table 1 we selected 10 bright (f60/i > 1 Jy) 'isolated' sources (Table 3) around 
which there are no other sources within twice the Gaussian radius. They therefore do 
not suffer the confusion problem seriously. The IRAS colors of these isolated sources will 
be studied in this section, which are presumably representative for the discrete sources 
in general. 

At the position of each isolated 60/zra source, we did source extractions on the other 
3 maps (i.e. the 12/zra map, 25 \im map, and the 100/zra map) using single elliptical 
Gaussian fittings. For each source and each map the source extraction is repeated 3 
times, with the size of the fitting area equal to 3' X 3', 4'. 5 X 4'. 5 and 6' X 6' respectively. 
All sources are detected in each of the three maps at least in two of the three runs with 
different fitting areas. In Table 3 we give their fluxes and the uncertainties in the IRAS 
bands, calculated using similar formulae to those given in Eq's (l)-(6). Some sources 
were detected only in two of the three Gaussian fittings (again with the size of the fitting 
area equal to 3' X 3', 4'. 5 X 4'. 5 and 6' X 6' respectively) at wavelengths other than 60/zra. 
In these cases, the corresponding flux and the uncertainty are calculated from the results 
in the two successful fittings, and are marked by a flag The uncertainties include 
the calibration uncertainties of the corresponding maps, which are at 11%, 9%, 4% and 
10% levels for the 12/zra, 25/zra, 60/zra and 100/zra maps respectively (Xu &: Helou 1994). 
The columns of Table 3 are arranged as follows: 
Col.(l): source ID in Table 1; 
Col. (2): Line 1: 12/zra flux /12/i in Jy; 

Line 2: uncertainty of /i2^; 
Col. (3): Line 1: 25/zra flux in Jy; 

Line 2: uncertainty of /25^; 
Col. (4): Line 1: 60/zra flux in Jy; 

Line 2: uncertainty oi feo^] 
Col. (5): Line 1: 100/zra flux /ioo/i in Jy; 



TABLE 3. 
IRAS Fluxes of Isolated Sources 



(1) 


(2) 


(3) 


(4) 




(5) 


ID 


fl2/i 


f25/i 


f60/i 




fl00/i 




unc 


unc 


unc 




unc 




Jy 


Jy 


Jy 




Jy 


1 


0.51 


0.47 


3.06 




13.43 




0.14 


0.07 


0.18 




1.62 


2 


0.11: 


0.21: 


2.08 




4.08 




0.01: 


0.04: 


0.09 




0.46 


3 


0.35 


0.44: 


1.53 




5.46 




0.04 


0.07: 


0.13 




0.70 


6 


0.21: 


0.36 


1.57 




4.43 




0.04: 


0.07 


0.08 




0.56 


10 


0.30 


0.43 


1.88 




9.08 




0.04 


0.07 


0.16 




1.48 


12 


1.12 


1.13 


5.46 




23.86 




0.24 


0.18 


0.54 




4.81 


27 


0.28 


0.40: 


2.58 




10.44 




0.08 


0.04: 


0.22 




1.55 


31 


0.30: 


0.82 


6.52 




18.71 




0.05: 


0.09 


0.34 




3.59 


38 


0.27 


0.34 


1.58 




5.22 




0.04 


0.05 


0.15 




0.98 


39 


0.32: 


0.49 


2.81 




10.84 




0.05: 


0.06 


0.15 




1.26 


Line 2: 


uncertainty of /ioo/i 










Fig. 6 is 


an IRAS color- color diagram (log (£12^/^25^) vs 


• log (f 60 


i/j/fioo/j)) of these 


isolated sources, which shows a 


clear anti- correlation 


between the colors. This anti- 


correlation h 


as been found for 


galaxies (Helou 1986) 


, for 


regions 


surrounding very 
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massive stars (Boulanger et al. 1988), and for a complete sample of small areas within 
the M31 disk (Xu and Helou 1994). The trend found by WS87 that regions of high 
60/zra to 100/zra intensity ratio correspond to relative depressions in the 12/zra to 100/zra 
intensity ratio may also be of the same nature (Helou, Ryter &: Soifer 1991). 

The FIR color ratio (f60/i/fioo/i) of the sources varies in the range from ~ 0.2, 
which is slightly higher than the color ratio of the diffuse emission (0.17, see Xu &: 
Helou 1994), to 0.5, with a mean of 0.30±0.04. The means of other two IRAS colors 

are < fi 2/i /W >= °- 14 ± °- 02 and < W/W >= °- 19 ± °- 02 - 

The mean < f60/i/fioo/i >= 0.30 ±0.04 corresponds to a color temperature of 26±1 
K [y 2 emissivity law). Cox &: Mezger (1988) and Rice et al. (1990) also found that for 
sources in the Galaxy and in M33 the average color temperature is about 30K, rather 
than 40K as assumed in WS87 and in other early papers (e.g. Cox et al. 1986). On 
the other hand, it appears that on average M31 sources are probably cooler than the 
M33 sources which have a mean f60/i/fioo/i = 0.40 ± 0.11 (Rice et al. 1990), although 
the difference is not significant because of the large error bars. This possible difference 
would be consistent with the optical studies of M31 HII regions which show: (1) the 
lack of bright regions and a possible excess of faint regions compared to other nearby 
spiral galaxies (Kennicutt et al. 1989; WB92), and (2) the giant HII regions in M31 
being more diffuse than the ones in the Galaxy and in M33 (Kennicutt 1984). It might 
also be due to the deficiency of the Very Small Grains in M31 disk (Xu &: Helou 1994). 

5. Luminosity of the HH-region-associated FIR component and the present- 
day star formation rate 

The FIR sources in M31 coincide generally with giant-HII or HII- complexes. This 
agrees fully with the results in the literature. Cox &: Mezger (1988) and Perault et al 
(1989) used the longitudinal profiles of Galactic emission in IRAS bands to separate 
the FIR sources and the diffuse emission in the Galaxy, and found that all the sources 
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are associated with bright Giant Molecular Clouds/HII complexes, spiral arm segments 
and the Galactic center. Rice et al. (1990) detected 19 sources in the IRAS maps of 
M33 (3' X 5' resolution), all coincident with HII regions cataloged at optical and radio 
wavelengths. For the Large Magellanic Cloud, with IRAS maps of resolution of ~ 3' X 3', 
Xu et al. (1992) found that the probability of coincidence between bright FIR sources 
and bright HII regions is better than 90 percent. On the other hand, it is very likely 
that some HII regions which are faint and do not cluster around any giant HII regions 
will not appear as bright FIR sources, namely that part of the HH-region-associated 
FIR emission may not be represented by the FIR sources. 

We try to estimate quantitatively how much the HH-region-associated component 
contributes to the total 60/zra emission of M31, through a comparison between the 60/zra 
sources and the HII regions detected in the CCD Ha observations of WB92 for some 
regions in the north-east part of M31. The idea is that such a comparison will yield 
a mean of the ieo^/Ha ratio for HII regions. Multiplying this ratio with the total Ha 
flux of M31 HII regions, an estimate of the total 60/zra flux of the HH-region-associated 
component can be made. Of the sources presented in Table 1, we find 12 cases which 
fall in the region covered by the Ha observations. For these 12 we show the sum of the 
Ha fluxes (uncorrected for internal extinction) of HII regions located within the region 
(defined by the major and minor axes and the P.A.) of the corresponding 60/zra source. 
Since the resolution of the Ha observations is much higher than that of our IRAS maps, 
an FIR source usually encloses many HII regions, of which we list only the brightest 
few. A weighted mean (weight=l/<r 2 ) of < f 60/i /f(Ha) >= 1.5(±0.2) 10 12 Jy/(erg cm" 2 
s _1 ) is found for these sources. 

It should be pointed out that some of the diffuse Ha emission which is not 
catalogued as extended structures in WB92, and thus is not included in the above 
calculation, so that the Ha fluxes may be underestimated and the ratio overestimated. 
According to Walterbos and Braun (1994), the diffuse emission contributes 40% of the 
total Ha emission of M31, some of which is included in the HII region catalogue as 
extended structures, and is not therefore missed by our estimate. On the other hand, 
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the Ha emission associated with the 60/zra sources is clearly dominated by a few very 
bright discrete HII regions. We therefore conclude that the error due to the omission 
of some diffuse Ha emission cannot be very significant. There are eight more 60/zra 
sources in Table 1 in the north-east half of M31, whose IDs are flagged with stars. As 
noted in Table 1, they are either completely or partially outside the Ha survey. 

Walterbos and Braun (1994), extrapolating from their survey over the NE part of 
M31, estimate that a total Ha luminosity of 1.2 10 40 erg s _1 (uncorrected for internal 
extinction) is radiated from the star formation regions in M31. This includes the flux 
from the diffuse ionized gas (DIZ) found near the regions of star formation which is likely 
to be photo-ionized by massive stars (Walterbos and Braun 1994), but excludes the 
nucleus and the diffuse component in the inner part of the M31 disk which appears not 
associated with star formation regions (Devereux et al. 1994). Devereux et al. (1994) 
find a higher Ha luminosity: excluding the nuclear emission and the diffuse emission in 
the inner disk, a total Ha luminosity of 1.8±0.5 10 40 erg s _1 (uncorrected for internal 
extinction) is observed from the 'star formation ring' in the M31 disk. However, this 
latter luminosity includes the emission from [Nil] lines. Using the spectral data of eight 
Sab - Sbc galaxies in Kennicutt (1992), we find a mean Ha/(Ha + [Nil]) flux ratio 
of 0.61±0.13. Hence the result of Devereux et al. (1994) implies a total 'pure' Ha 
luminosity of 1.1±3 10 40 erg s _1 radiated from M31 star formation regions, consistent 
with the estimate of Walterbos and Braun (1994). In the following we will use this value 
as the estimate of the total Ha flux from M31 star formation regions. 

At the adopted distance of 690 kpc, the Ha luminosity of 1.1 10 40 erg s _1 
corresponds to a flux of 1.9 10 -10 erg cm -2 s _1 . Multiplying this Ha flux with the 
mean < feo^/ f(Ha) > ratio found for the 60/zra sources, we estimate that a 60/zra flux 
of 289 Jy, with a relative uncertainty of ~ 40% (including both the error of the Ha 
luminosity and that of the < f60/i/f(Ha) > ratio), is due to dust associated with HII 
regions. This is 49±19 percent of the total 60/zra flux of M31, and twice the sum of 
f60/im of all sources listed in Table 1 and Table 2. Using this result and assuming that 
the mean IRAS colors of the isolated sources found in Section 4 represent those of the 
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HH-region-associated component in general, we estimate the fluxes of the HH-region- 
associated dust in the other three IRAS bands which, together with f60/i, are presented 
in Table 4. 



TABLE 4. 
IR Fluxes and Luminosities 
of HH-region-associated Component and of Diffuse Component 

/12/i /25/i /60/i /100/i fir" Lir 6 

Jy Jy Jy Jy 10" 8 erg cm" 2 s" 1 10 9 L Q 

HII 37±15 52±21 289±115 974±390 2.17±0.87 0.72±0.29 

diffuse 131±51 79±42 306±118 2099±502 3.64±0.96 1.70±0.48 

total 41 168±49 131±36 595±28 3073±316 5.81±0.40 2.42±0.56 
a Integrated FIR flux between 40 — 120/zra. 
6 Integrated IR luminosity between 8 — 1000/zra. 
c Including the nucleus. 

d IRAS fluxes are different from Xu Sz Helou (1994). See the erratum (Xu Sz Helou 
1995). 



From the IRAS fluxes listed in Table 4, and applying the formula of Helou et 
al. (1988), we find that the HII-region-associated component contributes 2.17(±0.87) 
10 -8 erg cm -2 s _1 of the integrated flux in the wavelength range 40-120/zra, namely 
37±15% of the total flux in this wavelength range radiated from M31 (Table 4). In 
order to estimate its contribution to the total dust emission (8 — 1000 /zra) we need 
to extrapolate both the HII-region-associated component and the diffuse component to 
longer wavelengths beyond 120/zra. Assuming that the very small grains are only half as 
abundant in M31 cirrus as they are in Galactic cirrus (Xu &: Helou 1994) and applying 
a model based on that of Desert et al. (1990) with the ISRF intensity equal to that 

14 



in the Solar Neighborhood (Fig. 7), we find that the total dust emission for the diffuse 
component is 1.13(±0.32) 10 -7 erg cm -2 s _1 , 3.1 times that of its integrated flux in the 
wavelength range 40-120/zra. The average spectrum of M31 sources is extended to the 
same wavelength range (Fig. 7) by using another model based on that of Desert et al. for 
IR emission near an 05 star with the dilution factor X05 = 0.003 (i.e. dust heated at 
18 pc away from the star), and with the assumption that both the very small grains and 
the PAHs are only half as abundant in dust associated with M31 HII regions as they 
are in Galactic cirrus (PAHs may be destroyed in HII regions, Boulanger et al. 1988). 
We then find that the total dust emission of the HH-region-associated component is 
4.8(±1.9) 10 _8 erg/cm 2 /sec, 2.2 times that of its integrated flux in the wavelength range 
40-120/zra. In terms of luminosity in the wavelength range of 8 — 1000/zra, we conclude 
that M31 radiates 2.42(±0.56) 1O 9 £ , of which 7.2(±2.9) lO 8 ^© (i.e. 30 ± 14%) is 
due to the HH-region-associated component (Table 4). 

It is interesting to note that the total IR emission of M31 found here is very close 
to the result of WS87 who found a total IR emission of 2.6 lO 9 i0 . WS87 used a 
simpler two-component model with the warm component being a modified black body 
(assuming a v 2 dust emissivity law) specified by a single temperature of 40K, and the 
cool component a fixed spectrum with the l60/i/Iioo/i = 0.16 (T<j = 21K). On the other 
hand, the total IR luminosity of the HH-region-associated component is a factor 3.6 
times the luminosity of the warm component of WS87, which is 2 lO 8 i0 . The low 
luminosity of the warm component of WS87 is likely to be a consequence of the high 
color temperature (40K corresponding to feo^f fioo^ — 1) adopted for the warm dust. 

From the total luminosity of the HH-region-associated component, we can estimate 
the present-day star formation rate of M31 using the formula given in Walterbos (1988): 
Rir = 1 X 10 _3 LiR(HII)/tiR, where Rir is the star formation rate in Mq /yr, Lir(HII) 
the total luminosity of the HH-region-associated component, and tnt the time scale over 
which massive stars contribute to the heating of HH-region-associated dust. Assuming 
t m = 2 x 10 6 yr (Th ronson &: Telesco 1986), we find a present-day star formation rate 
of Rir = 0.36(±0.14) M Q /yr, well in the range of 0.2 - 0.5 M Q /yr found by Walterbos 
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(1988) using several star formation indicators. This is about an order of magnitude lower 
than the present-day star formation rate of the Milky Way Galaxy (Mezger 1988). 

Our results show that, despite such a low present-day star formation rate, a 
substantial part (~30 — 40%) of the FIR emission of M31 is due to the young high 
mass ionizing stars also responsible for the HII regions. For the Milky Way Galaxy, the 
contribution from the HH-region-associated component to the total FIR emission is only 
20-30% (Cox & Mezger 1988; Perault et al. 1989; Bloemen et al. 1990). The relatively 
high contribution from the HII-region-associated component to the FIR emission of 
M31 is probably due to a depression of the diffuse FIR emission of interstellar dust not 
associated with HII regions, which in turn is caused by 1) a deficiency of the non-ionizing 
UV radiation which is responsible for most of the heating of diffuse dust in an average 
spiral galaxy (Xu 1990), given that the 2000A-to-blue flux ratio of M31 is one of the 
lowest in a large sample of nearby galaxies (Buat &: Xu 1995); 2) a deficiency of diffuse 
dust in M31, corresponding to a face-on V-band optical depth of Ty ~ 0.3 (Paper II) 
which is more than a factor of 2 lower than the average face-on optical depth of Sb/Sc 
galaxies (Xu &: Buat 1995). On the other hand, our results, showing that about 60% of 
the FIR (40 — 120/zra) emission of M31 is due to diffuse dust not associated with HII 
regions, do not support the suggestion of Devereux et al. (1994) that the FIR emission 
of M31 is predominantly due to heating by ionizing stars. 

Summary 

The high resolution (HiRes) IRAS maps of the Andromeda Galaxy (M31) are 
used to investigate the morphology of the infrared (IR) emission and the properties of 
discrete far-infrared (FIR) sources in its disk. Very thin and bright FIR arm segments 
are shown in these maps, which have similar structure as the HI gas but with much 
enhanced arm/inter-arm contrast, typically a factor of ~ 5 on the 60/zra image. We 
identify 39 unconfused sources (excluding the nucleus) and 14 confused sources (7 pairs) 
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at 60 /xm by direct Gaussian fitting to the image, with careful treatment of the source- 
confusion problem. The IRAS colors of 10 bright isolated 60/zra sources are studied, 
as a presumably representative sample of the discrete source population. They follow 
the well known anticorrelation between f60/i/fioo/i and fi2/i/f25/i, with mean IRAS colors 
< fw/Uon >= 0.14±0.02, < fasjx/feo,, >= 0.19±0.02, and < f 80 ^/fioo^ >= 0.30±0.04. 
All sources except the nucleus coincide with optical HII regions. A comparison with Ha 
observations shows that the total luminosity (in the wavelength range 8 — 1000/zra) 
associated with HII regions is 7.2 ±2.9 lO 8 i0 , namely 30 ±14% of the total IR emission 
of M31. This is a factor of 3.6 greater than the IR luminosity of the warm component of 
Walterbos and Schwering (1987). However, our result does not agree with the suggestion 
of Devereux et al. (1994) that the FIR emission of M31 is predominantly energized by 
high mass stars. About half of the IR luminosity of the HH-region-associated component 
emerges from discrete FIR sources. The luminosity translates into a present-day star 
formation rate of 0.36±0.14 Mq /yr, about an order of magnitude lower than that of 
the Milky Way Galaxy. 
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FIGURE CAPTIONS 

Figure 1. Gray-scale map (overlaid with contours) of M31 at 12/zra. The beam size is 
1'.7 (the same in the following maps). The contour levels are 0.3, 0.9, 2.7 MJy/sr. 

Figure 2. Gray-scale map (overlaid with contours) of M31 at 25/zra. The contour levels 
are 0.3, 0.9, 2.7 MJy/sr. 

Figure 3. Gray-scale map (overlaid with contours) of M31 at 60/zra. The contour levels 
are 0.3, 0.9, 2.7, 8.1 MJy/sr. The unconfused sources (Table 1) are represented by the 
white crosses, each showing the location, size, and the orientation of a single source. 
The open diamond represent the confused sources (Table 2), each standing for a single 
component in a pair in Table 2. 

Figure 4. Gray-scale map (overlaid with contours) of M31 at 100/zra. The contour 
levels are 1.2, 3.6, 10.8 MJy/sr. 

Figure 5. Comparison between the 60/zra surface brightness distribution (solid curve) 
and the HI column density distribution (dotted curve) along two cuts both passing 
through the center of M31. The HI column density is plotted with arbitrary units. 
Panel a: along the cut with the position angle P. A. =50°, in the NE half of M31. Panel 
b: along the same cut, in the SW half of M31. Panel c: along the cut with the position 
angle P.A.=24°, in the NE half of M31. Panel d: along the same cut, in the SW half of 
M31. 

Figure 6. #(60,100) = l60/i/Iioo/i v.s. R(12,25) = I^/hsu diagram of 10 isolated 
sources (see the text). 

Figure 7. The IR spectra of M31 sources and of the diffuse component. The solid line 
is the spectrum calculated using a model based on that of Desert et al. (1990) for the 
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IR emission near an 05 star with the dilution factor X Q 5 = 0.003, and assuming both 
Very Small Grains (VSGs) and PAHs are only half as abundant as they are in Galactic 
cirrus. The dashed line is to extrapolate the spectrum of the diffuse component to the 
whole IR wavelength range using a model based on that of Desert et al. (1990) for cirrus 
with the intensity of the IRSF equal to that in the Solar Neighborhood, and assuming 
the VSGs are only half as abundant as they are in Galactic cirrus. 
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